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Abstract 


The angular momentum independent statistical decay model is often 
applied using a Monte-Carlo simulation to describe the decay of prefrag- 
ment nuclei in heavy ion reactions. This paper presents an analytical 
approach to the decay problem of nuclei with mass numbers less than 60, 
which is important for galactic cosmic ray ( GCR ) studies. This decay 
problem of nuclei with mass number less than 60 incorporates well-known 
levels of the lightest nuclei (A < 11) to improve convergence and accu- 
racy. A sensitivity study of the model level density function is used to 
determine the impact on mass and charge distributions in nuclear frag- 
mentation. This angular momentum independent statistical decay model 
also describes the momentum and energy distribution of emitted particles 
(n, p, d, t, h, and a) from a prefragment nucleus. 


Introduction 

The description of heavy ion fragmentation is important for assessing damage to human cells and 
microelectronic equipment during spaceflight and on high-altitude airplane flights. Currently, few mea- 
surements of fragmentation cross sections are available because many collision pairs of interest exist for 
assessing the full galactic cosmic ray (GCR) spectrum interacting with various materials of interest. 
Other reasons include the unavailability of heavy ion accelerators and the high cost of measurement 
programs. Physical models (ref. 1 ) of heavy ion reactions must instead be used to provide databases for 
radiation transport codes (ref. 2) from which damage assessment may be made. 

Heavy ion fragmentation is usually described as a two-step process of abrasion-ablation. Abrasion 
describes the geometric overlap of two nuclei and the accompanying loss of material (nucleons) in the 
overlap region, with the size of this region dependent on the impact parameter of the collision. The large 
piece of nuclear matter remaining after the abrasion step is called the prefragment nucleus. Ablation 
occurs when the prefragment nucleus is in a state of excitation and decays by particle emission to a sta- 
ble configuration. The abrasion step is described by using geometric pictures (refs. 1 and 3) or optical 
models (refs. 4 and 5), and more recently with quantum models (refs. 6 and 7). The ablation step is 
described in the geometric picture (ref. 1) with the Rudstam charge distribution formulas and in the 
optical model approach with Monte-Carlo codes (refs. 8 and 9) for the statistical decay of an equili- 
brated system. 

Development of nuclear databases for GCR studies requires reaction models that are computation- 
ally efficient and accurate. References 6 and 10 detail efficient computational procedures for describing 
the abrasion step of the reaction. The present report describes the analytic solution to the statistical 
decay model (refs. 11—13) and a sensitivity study of available nuclear level density models. 

Statistical decay codes are often applied to a large number of nuclei, including ones with very large 
mass number A (A > 60), which are not of interest in GCR studies. This report focuses on a restricted 
range of nuclei typical of the prefragments formed during GCR fragmentation and uses analytic meth- 
ods that are easily coupled to abrasion and excitation energy models. This approach incorporates the 
well-known levels of the lighter mass nuclei (A < 12) and applies the statistical decay model to heavier 
mass nuclei. The statistical approach examines the decay (or evaporation cascade) of heavier mass pre- 
fragments (A > 1 1), along with the decay at higher energies of some lighter mass nuclei where the decay 
properties are not well-known. The energy and momentum spectra of the emitted particles ( n , p, d, t, h, 
and a) during decay are also represented consistently in the same statistical model. 

This report also describes the analytic solutions for fragment formation from an initial excited 
nucleus and the energy and momentum spectrum of the emitted or evaporated light ions. Physical inputs 
are then described, and numerical examples of the models are presented. 



Statistical Decay Model for Residual Nuclei Formation 

In the classical theory of statistical decay, the probability that a nucleus a with excitation energy E* a 
will decay by emitting a particle j with kinetic energy E ■ is given by (refs. 1 1 and 13) 
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where gj is the statistical weight (listed in table 1), c b j is the formation cross section, \ij b is the 
reduced mass M of j and the residual nucleus b in a — > b + j , and S“ is the separation energy given by 


S a j=M j + M b -M a (2) 

(We use subscripts j, k, l, ... to label the light ions emitted in a decay and superscripts a, b, c, . . . to 
label the parent and daughter nuclei.) In equation (1), w Q is the level density of the residual nucleus at 
the residual excitation energy E b = E a - S** - E .. Also, in equation (1) F ■ is the total probability for 
emitting j from a given by 

^ /j 

F a . = \ Ea ~ j P a AE) dE (3) 
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which satisfies the condition 

I F ° = 1 < 4 > 

j 

where the summation is over the types of particles which may be emitted at a given excitation energy. 
The actual values of input parameters for describing the probability function are described in the “For- 
mation Cross Sections and Decay Probability” section of this report. The present development is inde- 
pendent of these values; others could be used. Mass conservation and charge number conservation lead 
to 
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Z b + Zj 

(6) 


The probability of emitting j and forming a stable residual nucleus b depends on the residual exci- 
tation energy and the separation energies of b, denoted by S k . If the excitation energy is below the min- 
imum separation energy of the residual nucleus, no further decays are possible, which implies a 
stable b: 


E l = E a -S°j-Ej S min[S*] (7) 

For an unstable b, 

E* b = E* a -f j -E j > mm[Sj] (8) 
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For charged particle emission the conditions in equations (7) or (8) must be modified to include the 
Coulomb barrier. The total probability F a . at excitation energy E a is divided into two parts correspond- 
ing to formation of a stable or unstable residual nucleus as 


F a j = G a j+H a j 


(9) 


where the probability of forming a stable residual nucleus is given by 
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and the probability of forming an unstable residual nucleus is given by 



mints*] a 
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For sufficiently high initial excitation energy E a several to many decays may occur; this occur- 
rence leads to an evaporation cascade. In the statistical model the individual steps are treated as largely 
independent and connected only through the conservation of energy, mass, and charge. The evaporation 
cascade can be represented by using a master equation (ref. 12), with the development herein represent- 
ing a pertubative solution to this equation. 

* b 

The possibility of a second decay occurs for E b > min [S^J , and the probability of forming a stable 
nucleus from b is given by 
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with similar expressions for higher order terms. We approximate the second and higher order terms in 
the evaporation cascade by introducing the average energies of Ej in emission to form stable residuals 
given by 



f f l b EP“(E ) dE 
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and unstable residuals given by 
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The probability of forming a stable residual nucleus with A c = A a - Aj~ A k and Z c = Z a ~Z .-Z k 
after two decays is then approximated by 



nh PL * /) ^ n — H 

Gjt~H a j(E a ) G k (E a -Sj-Ej ) 

(15) 

and a stable residual nucleus with A, = A - A A-- A, and Z, = Z -Z-Z.-Z, 

after three 

decays by 
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Gjki = Hj(E a ) H b k (E a - Sj - Ej ) G c ,(E a - S . - S k - Ej - E t ) 

(16) 


with similar expressions for higher order terms. The number of terms in the cascade series depends on 
the initial excitation energy and the number of decay particles j allows. Early Monte-Carlo studies 
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(ref. 13) suggest that for nuclei in the mass range of the GCR (A < 60) the probability of evaporation of 
particles heavier than 4 He is small. Nuclear structure effects on the separation energies reduce the likeli- 
hood of emission of the mass 2 and 3 ions, especially for A > 20, although structure effects may play a 
role in some cases. Table 2 presents the number of terms in the decay series for emission of n, p, and a 
with emission of n, p, d,t,h, and a, which is clearly a slowly converging series for large prefragment 
mass (A > 20) and high initial excitation energy. 

Spectrum of Light Ions 

The spectrum of emitted ions j is found by summing the contributions of the j ions from single, 
double, and so forth, evaporations until the initial excitation energy is exhausted. Keeping track of the 
excitation energy available in the decay probability functions by introducing the notation 

P a j(E* a >Ej ) = P a j(Ej) (17) 

for equation (1) allows the cumulative spectrum to be written as 


dN ; 
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The higher order terms in equation (18) correspond to high initial excitation energy. At high excitation 
energy, the statistical decay properties (as discussed in the section “Results and Discussion”) will 
change slowly. An approximation of equation (18) similar to the one used in equations (15) and (16) is 
then found as 


dNj 

~dE 


= P a j(E’ a , E) + X H“(EJ p)(E a - S° k - E k , E) 
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H -H 


+ ^T H k( E a) H ,( E a- s k - E k) Pj(.E a -Sl-Sl-E, -£,,£) + 


k l 


(19) 


Assuming isotropic emission, the momentum distribution in the rest frame of the initial nucleus is found 
as 


dN . i i dN . 

I = — — — — l (20) 

dp 4n pE dE 

In heavy ion fragmentation the prefragment nuclei are produced with a velocity or momentum recoil 
spectrum in the first stage of the reaction. The recoil spectrum is expected to be well represented by a 
Gaussian function and is denoted by /(v - v Q ), where v Q is the velocity vector of the primary nucleus, 
Vq = 0 for target fragmentation, and v Q ~ \ p , where \ p is the projectile velocity vector for projectile 
fragmentation. The invariant momentum distribution in the lab or projectile rest frame is then 

, dN. r dN. 

£ ^=i rfv/(v - v « )£ ^ (21) 
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where primed values refer to the laboratory or projectile frame and 

Ej = Y v (Ej - |3 V p'j cosGy) (22) 

where y v and (3 y are functions of v and v () . 

Formation Cross Sections and Decay Probability 

Following Dostrovsky (refs. 13 and 14), the energy dependence of the formation cross sections, the 
effects of tunneling, and the Coulomb barrier are simulated by using the following method: 

For neutrons, 


= nR j a n 
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with 
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and 


p. 


— (1.22 A ^ -0.05) 
ct 0 


For charged ions, 


R. = 1.5 A 


- 1/3 
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with the Coulomb barrier given by 


v j b = 


Z b Z j e 
R b + Pj 


(25) 


(26) 


(27) 


(28) 


The parameters a ., k j, and p 7 are given by Dostrovsky (ref. 14) for Z b > 10 and are extrapolated lin- 
early to lower Z b herein. 

In the statistical model, the level density at E - S . - E . is approximated by (ref. 13) 

w J J 


b * a b * -(■ Ei + SA/ \ 

w b 0 (E -S a j-Ej)~w b 0 (E )e J J J 


where t . is the nuclear temperature defined by 
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(29) 


(30) 


The decay probability function is then rewritten in the form 


a a a b * -(£.• + Sy)/T, 

P a j(E.) = A a j(E r B^)w b 0 (E a )e J J J Q(Ej- K.V hj ) 


(31) 
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where 
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and 
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The functions Fj, G“, and H 6 * then have solutions of the form 


a a b * r E 2 a -{E. + fy/T* , 


(32) 


(33) 
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with E l = B a j, if E { < b“ for j & n. Similarly, the average energy of the emitted particles Ej and Ej 
can also be found in closed form. 


In calculations, the separation energies and masses of the particles are evaluated by using experi- 
mental values from reference 15. A large number of decays may occur for high excitation energies. The 
function Gj approaches zero at high E * , whereas H j becomes almost constant. Calculations are shown 
in the section “Results and Discussion;” however, first we discuss models of the level density function. 


Models for Level Density Function 

The level density function can be derived from an equidistant spacing model that assumes all one- 
particle states are equally spaced. The level density then takes the form (ref. 1 6) 


w Q (A, Z, E ) 




* 5/4 

12 E a 


1/4 


where a is the level density parameter and S the entropy, which is expressed as 


(36) 


S = 


2 *JaE 


(37) 


In the Fermi gas model, the level density parameter is related to the spacing parameter d by (ref. 16) 


a = 


n 
6 d 


with d expressed in terms of the Fermi energy s F and mass number A by 

2£rr 


d = 


3 A 


(38) 


(39) 
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The values used for the level density parameter found in the literature often range from a = A / 8 MeV 1 
to a = A/20 MeV -1 . 

To account for important deviations from the equidistant spacing model such as pairing effects, 
shell structure, and energy gaps, the excitation energy in equation (36) is often replaced by an effective 
excitation energy, and the level density parameter is allowed to have a dependence on the excitation 
energy. The model of Ignatyuk, Smirenkin, and Tishin (ref. 17) uses the form 


5 2 = 


4a 


E -A p -A, 


+ A£ 


shell 



(40) 


where A p and A n are the proton and neutron pairing corrections, respectively, and Ais shell is the 
ground state shell correction given by 


^ shell ~ M cxp M LD 


(41) 


where M exp is the experimental mass and M LD is the mass corresponding to the liquid drop model. 
The pairing correction is taken as 


n or p 


12 MeV 

• A 1/2 

. 0 


(if n or p is odd) 
(if n or p is even) 


(42) 


The level density parameter is written as 

a = A(cj + c 2 A) 


(43) 


Several modifications to equation (40) have also been considered (ref. 16), including studies where the 
parameters c,, c 2 , and y are fit to slow neutron and proton resonances for most nuclei or a restricted 
range of mass numbers. The equation of state corresponding to equation (40) is 


* 



1 +A£ 


1-^M 2 


shell 


(44) 


Several other corrections to the level density model and decay parameters can also be considered. 
These include use of a limiting temperature to restrict states with short-lived lifetimes that are not 
expected to decay statistically. (See ref. 18.) Also, temperature-dependent masses and Coulomb barriers 
are often introduced at high excitation energies. (See refs. 13 and 19.) These corrections can be studied 
in this model by modification of input parameters. 


Results and Discussion 

The use of an average residual energy of the daughter nuclei to replace the integral over the energy 
spectrum of emitted ions is the main approximation used to reduce the evaporation cascade to a more 
convergent summation of terms. To test this approximation, figure 1 presents calculations using the 
approximate expressions and approximations using numerical integration for the evaporation spectrum 
of n, p,d , and a ions from 19 F. Terms corresponding to three decays at several initial excitation ener- 
gies are shown. The approximate solutions work well, and similar results are found with other nuclei. 

To increase the accuracy and convergence of the evaporation cascade, the functions Gy, and 
H U j for lighter mass nuclei (A < 11) are replaced by experimental values (refs. 20-22) up to the highest 
known excitation energies where the statistical model is used at higher excitation energies. The statisti- 
cal model agrees reasonably well at low excitation energy for higher mass numbers. Figures 2-5 show 
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the functions G and H versus excitation energy for a representative number of nuclei. In most cases the 
functions G a are small above E = 50 MeV, and the functions Hj are fairly constant above 100 MeV. 
This behavior allows the use of G°j ~ 0 in the evaporation cascade at high excitation energies 
(>100 MeV) and improves convergence in the code. 

To consider the effect of the level density model employed when calculating predicted fragmenta- 
tion cross sections, the abrasion-ablation model is used to consider the effects on elemental and mass 
distributions for realistic cases. The abrasion cross sections are calculated in the abrasion model of 
Townsend et al. (refs. 5 and 23) with numerical techniques described by Cucinotta, Townsend, and 
Wilson (ref. 10). Here the elemental distribution is given by 

°z F = X £ £ °v, p* °abr(V’ V’ ■*>*) . < 45) 

a f V v 

and the mass distribution by 

= X X X a F F * °abr(V' V’ v) (46) 

Z F V V ’ 

where the parameter a F F * represents a cumulative G a ^ function that results in the formation of a 
fragment F from the prefragment F . 

In the present calculations, the excitation energies of the prefragments are from the NUCFRG code 
(refs. 1 and 24), which gives two excitation energies for each prefragment corresponding to surface and 
frictional excitation energies at two trajectories. Figures 6-10 show elemental and mass distributions for 
several projectile nuclei interacting with 12 C targets. References 25-27 present experimental data for 
elemental production cross sections. Calculations are shown that use the level density function of 
equation (36) in the equidistant spacing model with a = A/8 MeV -1 . Calculations including pairing 
and shell corrections using the effective excitation energy (eq. (43)) noted in reference 17 are also 
shown. The equidistant spacing model without the level density corrections is in poor agreement with 
the experiment for most odd Z fragments. In contrast, use of the effective excitation energy leads to 
fairly good agreement for the projectiles studied for both even Z and odd Z fragments. Large differences 
between the two models for the level density are also seen in the mass distribution of fragments. 

The present abrasion model needs further corrections for elastic fragmentation and excitation of low 
level states which decay through gamma-ray emission to improve the agreement for fragments 
produced with masses near that of the projectile. Here the NUCFRG model does not consider the actual 
separation energies, and with removal of a few nucleons the excitation energies may be mismatched 
with realistic models of level spectrum. Microscopic approaches to define the excitation energy spec- 
trum of the prefragment nuclei must be developed to improve predictive capability and identify short- 
comings in the statistical decay model. 

The calculations in this report for fragmentation of projectile nuclei with A P < 28 have a CPU time 
under 5 minutes on a VAX-4000 series computer. As mass removal increases, prefragment excitation 
energies exceed 200 MeV and the number of prefragments grows dramatically. The required CPU time 
for 40 Ar is about 2 hours and for 56 Fe about 8 hours. Here a Monte-Carlo simulation of the evaporation 
cascade would be more useful. However, the convergence of the analytic approach may be improved by 
neglecting d, t, and h emissions, which greatly reduces the number of cascade terms for large numbers 
of decays (high excitation energies). The importance of d, t, and h emissions is expected to be reduced 
for A f * > 20 because of the increased separation energies of these ions compared to p, n, and a separa- 
tion energies. Figure 10 presents elemental distributions for three or six decay particles for 20 Ne, 40 Ar, 
and 56 Fe fragmentation. The truncation to three decay particles (n, p, and a) appears to be fairly 
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accurate for heavier ions such as 40 Ar and 56 Fe. A hybrid model using six ions for excitation energies 
below about 100 MeV and three ions for excitation energies above 100 MeV for the heavier prefrag- 
ments is clearly suggested. 

Concluding Remarks 

Improvement in the database for heavy ion fragmentation cross sections is the primary area of phys- 
ics research for development of radiation transport codes for cosmic ray studies. The correct physical 
description for developing models of these reactions is the two-step model of abrasion-ablation. Past 
research has not focused on describing the ablation step when developing heavy ion fragmentation data 
bases. This report shows that an analytic formalism of the statistical decay model provides reliable 
results when using the abrasion- ablation model. Furthermore, the model shows that the nuclear level 
density used in the calculations has a major effect in predicting secondary mass and charge production 
cross sections. The approach developed in this report can be extended to aid the development of a com- 
plete database for galactic cosmic ray radiation transport codes. Improving physical inputs to the calcu- 
lations should be considered, and improved convergence methods for the evaporation cascade at large 
excitation energies in the mass 25-45 region should be developed. 


NASA Langley Research Center 
Hampton, VA 2368 1 -0001 
July 16, 1996 
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Table 1 . Light Ion Statistical Weights and Masses 


Ion j 

g j 

Mj , MeV 

n 

2 

939.55 

P 

2 

938.26 

d 

6 

1875.58 

t 

6 

2808.87 

h 

6 

2808.34 

a 

4 

3727.32 


Table 2. Number of Terms (tj or t 2 ) for Number of Decays With 3 or 6 Decay Particles 


Number of 
decays, N 

t j (3 particles) 

*2 (6 particles) 


1 

3 

6 

2 

2 

9 

36 

4 

3 

27 

216 

8 

4 

81 

1296 

16 

5 

243 

7 776 

32 

6 

729 

46656 

64 

7 

2187 

279 936 

128 

8 

6561 

1 679 626 

256 

9 

19 683 

10077 856 

512 
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(a) E* = 25 MeV. 

Figure 1. Energy spectrum of evaporated particles from 19 F at different levels of MeV excitation 
energy using numerical integration of equation (18) and approximate solution of equation (19). 
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(b) 25 Mg. 

Cl Cl 0^1 oc 

Figure 5. Calculations of Gj and H . versus excitation energy for Al, Mg, and Na. 
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(b) Mass production. 

Figure 6. Excitation energy cross section profiles for 20 Ne + 12 C at 0.6 GeV/amu. 
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(a) Elemental production. 



(b) Mass production. 

Figure 7. Excitation energy cross section profiles for 24 Mg + 12 C at 0.74 GeV/amu. 
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(a) Elemental production. 



(b) Mass production. 

ryq 1 r\ 

Figure 8. Excitation energy cross section profiles for 'A1 + C at 0.58 GeV/amu. 



(b) Mass production. 
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Figure 9. Excitation energy cross section profiles for J o + C at 3.65 GeV/amu. 
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(c) Cross sections for 56 Fe + 12 C at 0.52 GeV/amu. 


Figure 10. Concluded. 
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